Optimal reference states for maximum accessible entanglement under the local 

particle number superselection rule 



o 
o 

(N 






> 

o 
p 

o 

oo 

o 



X 



G.A. White, J. A. Vaccaro, and H.M. Wiseman 

Centre for Quantum Computer Technology, Centre for Quantum Dynamics, 

School of Science, Griffith University, Brisbane, 4 m Australia 

(Dated: January 19, 2009) 

Global conservation laws imply superselection rules (SSR) which restrict the operations that are 
possible on any given state. Imposing the additional constraint of local operations and classical 
communication (LOCC) forbids the transfer of quantum systems between spatially separated sites. 
In the case of particle conservation this imposes a SSR for local particle number. That is, the 
coherences between subspaces of fixed particle number at each site are not accessible and any 
state is therefore equivalent to its projection onto these subspaces. The accessible entanglement 
under the SSR is less than (or equal to) that available in the absence of the SSR. An ancilla can 
be used as a reference system to increase the amount of accessible entanglement. We examine the 
relationship between local particle number uncertainty and the accessible entanglement and consider 
the optimal reference states for recovering entanglement from certain systems. In particular we 
derive the optimal ancilla state for extracting entanglement for a single shared particle and make 
steps towards the optimum for general systems. We also show that a reference for phase angle is 
fundamentally different to a reference for the SSR associated with particle conservation. 

PACS numbers: 03.67.-a, 03.67.Mn, ll.30.-j, 03.65.Ta 



I. INTRODUCTION 

Some five decades ago Wick, Wightman and Wigner 
showed that conservation laws induce so called "supers- 
election rules" (SSR) that forbid the observation of co- 
herences between different eigenstates of the conserved 
quantity by any physical measurement [l| . A decade later 
Aharanov and Susskind demonstrated that this restric- 
tion can be alleviated with an additional system that acts 
as a reference for the coherences [2j. Such reference an- 
cillae are of particular interest when the dimension of the 
Hilbert space it occupies is not arbitrarily large compared 
to that of the system. In such cases the quantum nature 
of the reference system is important. We refer to such 
a reference system as a finite reference because of the 
finite dimension of the reference system's Hilbert space 
[36|. This is the regime that we are mainly concerned 
with here. 

There has been some laxity in the use of terms to dis- 
tinguish the classical and quantum regimes in the litera- 
ture. To avoid any possible confusion we shall reserve the 
term reference frame for a set of independent variables 
which, in principle, can be defined without an explicit 
relation to any physical object. The inertial reference 
frames in special relativity are an example. In contrast, 
we use the terms reference system and reference ancilla 
to describe a physical system whose observables are used 
as a reference for other physical systems. This distinc- 
tion between frames and systems is crucial in theories 
such as quantum gravity where external reference frames 
are artificial and a description of a physical system can 
only be made relative to other physical systems [3, |j| . 

The lack of the ability to observe coherences under a 
SSR implies that there are operations whose effect on the 
system are not physically detectable. Conversely, the ex- 



istence of such operations implies that a SSR operates. 
This converse situation is true even in the absence of a 
conservation law. For example, it can be shown that the 
set of non-detectable unitary operations form a symme- 
try group, G 0. In other words, G is the symmetry 
group which describes the invariance of the system un- 
der the SSR. Moreover, this invariance implies that the 
set of physically allowed operations {C(-)} is constrained 
to commute with the symmetry group G [6|, that is, 

0[fig)pf^ig)]^fig)0[p]f^ig) ^g e G , (1.1) 

for every state p of the system. Here T{g) is the unitary 
representation of the transformation g G G. In general, 
O represents an operation on an open system and the 
openness implies that any global conservation law may 
not hold for the system itself. This shows that SSRs can 
exist independently of conservation laws at least for open 
systems. 

In accord with previous work, we label the SSR associ- 
ated with the symmetry group G as the G-SSR. Alterna- 
tively, if G is characterized by a generator, we sometimes 
refer to the SSR by the name of the generator. For ex- 
ample, a unitary representation of U(l) for a system of 
identical particles is given by {e*^^ : < 9 < 27r} where 
the group generator N is the particle number operator, 
and so here the U(1)-SSR is equivalent to the particle 
number SSR. 

The situation is further enriched if, in addition to a 
SSR, we also impose the constraint of allowing only lo- 
cal operations and classical communication (LOCC) for 
a bipartite system. The physically allowed operations O 
must satisfy both Eq. (|l.ip and LOCC. Consider the bi- 
partite case with subsystems at two spatially separated 
sites A and B. Any local operation O = Oa <8) Is that 



satisfies Eq. (|l.ip has the property that 

{OA®tB)[f{9)pf\g)]=f{g){OA®^Bmf\g) (1.2) 

for all g ^ G, where O^ is a local operator acting at site 
/i and 1 is the identity operator. If the group has local 
representations {T^((7) : g g G} for each site /i where 
f{g)^fA{g)^fBig) then 
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For the case of a Lie group, this requires its local genera- 
tors to be additive. Thus the combination of LOCC and 
a SSR can induce a local SSR acting at each site. Such 
local SSRs are not necessarily fundamental in the sense 
of arising from a conservation law like that of charge, but 
are imposed by our interest in examining the effect of the 
LOCC constraint. Essentially, a local SSR arises when 
two sites lack a shared reference to break the symmetry 
represented by G at each site [a, S 0, B- However, the 
absence of the shared reference here is due to the LOCC 
constraint rather than any fundamental constraint. For 
example, the local U(1)-SSR due to the absence of a 
shared optical phase reference is imposed by the prohi- 
bition on transporting quantum fields between two sites 
under the LOCC constraint rather than the conservation 
of photon number. 

Interest in superselection rules and their associated ref- 
erence systems has revived in recent years particularly 
in relation to quantum information (9|. Bartlett et al. 
have shown that in the absence of a reference for spatial 
orientation the communication of classical and quantum 
information using spin particles is still possible provided 
the spins are entangled and the information is encoded 
in SSR- invariant subsystems [lOj. They have also shown 
how a private shared reference for spatial orientation can 
be used for secret communications |ll|. The effects of 



finite references has been studied in a variety of ways. 

For example, Bartlett et al. 10], Bagan et al. [lj| 
and Lindner et al. [13i] have examined protocols for es- 
timating the relative angle between two directions which 
are defined by finite dimensional spin systems. Further, 
the effect on quantum operations due to finite refer- 
ences was studied initially by van Enk and Kimble [ij| 
and shortly later by Gea-Banacloche [1^1 and, although 
specific details of this initial work attracted some crit- 
icism from Itano [ig] and Nha and Carmichael [l7|, it 
has been extended and generalized to quantum measure- 
ments, uncertainty relations and simultaneous measure- 
ments mainly by Ozawa ^, IH [13, Hll ■ 

Entanglement is also affected by the presence of a SSR. 
A SSR has an effective decohering effect that constrains 
the amount of entanglement in a system that is accessi- 
ble. A number of different terms have been used to la- 
bel the entanglement which is constrained in this way. 
Bartlett and Wiseman [6|], as well as others (see e.g. 
d:. 3, [23, 1231), refer to it as the accessible entanglement 
37| ■ For the special case of indistinguishable particles in 
the absence of a shared phase reference, we previously 



referred to the SSR-constrained entanglement as parti- 
cle entanglement in Refs. 0, [M] to distinguish it from the 
entanglement of the spatial modes occupied by the parti- 
cles. Our operational definition of particle entanglement 
makes use of a set of local quantum ancillary systems (or 
registers) which are not subject to the SSR. The accessi- 
ble entanglement in a system of shared identical particles 
is given by the maximum amount of entanglement that 
one can transfer from the system to the local ancillas by 
G-invariant operations that obey LOCC. It is important 
to note that the transfer operation involves only a single 
copy of the system state. However, once the entangle- 
ment resides in the ancillae it can be manipulated in the 
usual way free of the SSR. Further details of this non- 
asymptotic interpretation of particle entanglement can 
be found in the Introduction of Ref. '5| . 

In this paper we focus on the problem of finding the op- 
timal reference state that maximizes the accessible entan- 
glement of a system in the presence of the local particle 
number SSR. Under such a SSR, the pure state describ- 
ing a single particle which shared coherently between two 
sites is physically equivalent to an equal mixture of states 
representing a single particle at one site and a single par- 
ticle at the other. In general, the SSR constrains the 
accessible entanglement due to the unobservabilty of the 
coherences between eigenkets of different local particle 
number. The accessible entanglement in this case is the 
weighted average of the entanglement found in each sub- 
space of fixed local particle number. The entanglement 
lost because of the SSR can be recovered with a reference 
ancilla and it is this problem that we focus on here. In 
particular, we examine cases where the number of par- 
ticles in the reference ancilla is not arbitrarily large and 
we explore how effectively the reference ancilla increases 
the accessible entanglement in a variety of situations. In 
Section |TT] we calculate the optimal reference state for 
maximizing the accessible entanglement from a system 
containing a single shared particle where the number of 
particles in the reference system is fixed. We compare the 
amount of entanglement made accessible by the optimal 
state with that due to various other states in section IIIII 
and make steps towards finding optimal reference state 
for a general state of the system in section IIVI 



II. 



OPTIMAL REFERENCE STATE FOR A 
SINGLE SHARED PARTICLE 

A. Particle entanglement 



We begin by briefly recalling the definition of the parti- 
cle entanglement (or equivalently, the accessible entangle- 
ment) as defined by Wiseman and Vaccaro [3] . As shown 
above, SSRs can be induced by the lack of a shared refer- 
ence or a conservation law. The origin of the SSR is unim- 
portant for the main results of this paper. However, as 
the conservation of particle number was used in Ref. [7| , 
we assume the conservation also holds here. The obser- 



vation of coherences between states of different particle 
number requires operations that do not conserve particle 
number. The conservation of particle number therefore 
implies that phase coherences between subspaces of the 
Hilbcrt space corresponding to different numbers of par- 
ticles are unobservable. This means that phase shifts 
generated by exp(—iN9), where N is the particle num- 
ber operator and 9 a phase angle, are not detectable. 
The group of undetectable transformations is therefore 
the group U(l)= {exp{-iNe) : < < 2tt} oi phase 
shift operators, which induces the particle number super- 
selection rule U-SSR. In the case of bipartite systems, 
we consider the SSR which is induced by imposing the 
conservation of particle number at each spatial site. The 
corresponding group is Uab(1) = Ua(1) <8i Ub(1) where 
U^(l) is the group of operators that generate phase shifts 
at site fj,. We call this the local particle number superse- 
lection rule, i.e. the local U-SSR. 

Consider the pure state |^*-^-')ab of a system compris- 
ing a fixed number of indistinguishable particles shared 
between two spatially separated sites labeled A and B. 
We have included a superscript in the label of the state 
to indicate the total (fixed) number of particles the state 
represents. The effective state under the local particle 
number SSR is not |\E'^^')ab but rather a mixed state 
given by 
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where il^g " is a operator which projects onto the 
subspace representing n particles at site A and (N—n) at 
B. The entanglement of p^^-' is, by definition, the parti- 
cle entanglement of |*''^')ab- We note that the effective 
state, p^^\ is a mixture of a set of mutually orthogonal, 
entangled pure states n)(*B ~" |\E'^^^)ab- Its entangle- 
ment is thus the average of the entanglement of each 
member of the set. This means that the particle entan- 
glement is given by Q 
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Here S{g) is the von Neumann entropy of g, |^„ ■^)ab is 
the state of the system after detecting n particles at site 
A, and Em{\iP)) is the entropy of entanglement of the 



state \ip). In the current context EmU'iP)) corresponds to 
the entanglement of the spatial modes that are occupied 
by the particles. 



B. Particle entanglement of system and reference 
ancilla 

Two of us previously showed [3] that two copies of a 
shared single-particle system contained particle entangle- 
ment whereas (as noted above) each of the systems does 
not. Evidently one system behaves as a reference for the 
other. We further developed developed this concept in 
Ref. [a| where we showed that a shared particle ancilla 
can be used as a reference to increase the accessible the 
entanglement in another system. We now examine the 
problem of finding the state of the reference ancilla that 
yields the maximum value of the particle entanglement 
of the combined system and reference ancilla for a given 
state of the system. We confine our analysis to situa- 
tions where the system and reference ancilla are in pure 
states to avoid unnecessary detail. Also in accord with 
our assumption that the particle number is conserved, 
we constrain the total number of particles in each of the 
system and the reference ancilla to be fixed [3a] ■ Our fo- 
cus will initially be on the simplest case, that of a single 
shared particle in the state: 
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where the ket \n, to)ab represents the occupation of spa- 
tial modes at the sites A and B by n and m identical 
particles, respectively, in second quantization notation. 
It is straightforward to show that this state possesses no 
particle entanglement, i.e. £;p(|*(i))ab) = U,\^. Let 
the reference ancilla consist of M particles of the same 
type as the system and be prepared in the pure state 



M 
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n=0 



for which the coefficients c„ satisfy the normalization 
condition 



M 



{<^iM)^^iM}^^J2\Cn\^ = l 



(2.7) 



n=0 



but are otherwise undetermined. This ancilla state has 
been intentionally constructed to have no particle entan- 
glement itself, i.e. Ep{\^^^'^'>)ab) = 0. Any particle en- 
tanglement in the combined system and ancilla state, i.e. 
in 



|C(*^+l))AB = |*W)AB®|<i'(*')) 



AB 
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will therefore be due to the ancilla's ability to ameliorate 
the effects of the local U-SSR. 



Expanding |C^*'^"'""'^^)ab in terms of states with fixed 
local particle number yields 



and so, from Eq. (|2.3p . the particle entanglement is the 
average 



M+l 
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is a state of representing a total of n particles at site 
A and (M + 1 — n) at site B. For convenience, we omit 
the site labels AB here and below when they are implied 
by the context. The projection operator II 



(ri,Af+l-n) 
AB 
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of the same basic form as the one in Eq. (|2.4p but here 
it projects onto the subspace of states of the combined 
system plus ancilla with a total of exactly n particles at 
A and (A-f + 1 — n) at B. That is, it distinguishes between 
the sites but it is insensitive to the individual modes at 
each site. Making use of Eq. (|2.5p we find that 



|C(^^+i)) = ^(c„_i|l,0) ® In- l,M-n+ 1) 

(2.11) 
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for n = 0, 1, . . . , il/ + 1. To make expressions more com- 
pact here and in the following, we have introduced two 
extra coefficients c_i and cm+i whose values are zero, 
i.e. 



|c-iP = |cA/+i|2=0 



(2.12) 



We shall refer to Eq. (I2.12|) as the "boundary con- 
ditions". As before, the value of p„ ensures that 
(Cn \Cn ) = 1 and is easily found to be given 

by 



Pn = ■:^{\c7i-l\ + \Cn\ ) 
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for < n < M + l. Taking the partial trace of 
\C'n )(Cn -^ I over a basis representing states at site 
B yields 
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The entanglement of modes i^Mdw )) is given by 
the von Neumann entropy of Eq. (|2.14|) . i.e. 



EM{\Ci^'^'^)) = --^|c„_ip log,(^^) 
2p„ 2p„ 
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C. Conditions for the optimum reference 

The maximization of the particle entanglement in 
Eq. (|2.16p over the coefficients c„ is subject to the nor- 
malization of the ancilla state |<I>(*^)) given in Eq. (|2.7p . 
We note that the optimization can be performed with 
respect to |c„p rather than c„ since the particle entan- 
glement Eq. (|2.16p and the constraint Eq. (|2.7p are both 
functions of |c„p only. Let 

M 

F = (21og2 e)i?p(|C(*^+i)» - a(E W 1) (2.17) 

n=0 

be the auxiliary function where a is a Lagrange multi- 
plier. 

The coefficients of the optimal reference state are given 
by the extremum of Eq. (|2.17p and so satisfy 

BF 

g^ = -2[ln(|c„n + 1] + [ln(|c„_ip + |c„n + 1] 



= ln 
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+ [ln(|c„p + |c„+i|2) + l]-a 
(|c„-ip + |c„p)(|c„|2 + |c„+i|2) 
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for n — 0, 1, ... , M . This equation can be expressed in 
terms of a recurrence relation as follows 

/3|C„|4 = |C„|4 + |c„_i|2|c„+i|2 -I- |c„P(|c„_iP + |c„+in 

(2.19) 
for n = 0,1,..., M where we have set /3 = e". The 
boundary conditions in Eq. (|2.12p ensure the recurrence 
relation has the correct form for n = and n ~ M . 

To find the solution it is convenient to first rearrange 
the recurrence relation as 



C^n-l-ll 



(/3-l)|c„|4-|c„p|c„_i|- 

|c„|2-H|c„_i|2 



(2.20) 



Iterations of this starting from ri = lead to explicit ex- 
pressions for the coefficients in which |c„p is a polynomial 
in /3 of order n. The details are given in the Appendix. 
Unfortunately, the derivation of a solution from these 
expressions that satisfies the upper boundary condition 



0.07 





FIG. 1: Optimal reference state for a system comprising a 
single shared particle. The probability |c„p is plotted as a 
function of the number of particles n at site A for the nu- 
merical solution (crosses) and our ansatz (solid curve) for an 
ancilla with a total of M = 29. For the ansatz we used the 
values A — 1 and e = | which are correct to 0{M~^). 



|cA/+iP — and the normalization condition Eq. ()2.7|) 
does not appear tractable analytically. Nevertheless we 
note that the system state Eq. (|2.5p is symmetric with re- 
spect to the interchange of labels, \n,m)j^B '^^ |"^7?^)abi 
at each site A and B, which implies that if there is a 
unique optimal reference state, it will also be symmetric 
under the same operation. Indeed we also show in the 
Appendix using the polynomial expressions that 



|ca/- 



(2.21) 



for n = 0; 1, . . . , M. This result will be useful later. 

The optimal reference state for any given value of M 
can be determined numerically by computing the values 
of /3 and the set of coefficients {|c„p} that satisfy the re- 
currence relation Eq. (j2.20p . the normalization condition 
Eq. ([27ll and the boundary conditions Eq. (|2?T2l) . The 
approach we adopted was as follows: (i) make an initial 
guess of the value of /3, (ii) use the lower boundary con- 
dition |c_ip — and the recurrence relation Eq. (|2.19|) 
to compute the values of the variables |c^p = |c„p/|coP 
for n = 0, 1, . . . , M, (iii) use the computed values of |c^p 
and the normalization condition Eq. (|2.7p to determine 
the values of the normalized coefficients |c„p. This yields 



a set of normahzed coefficients |co| , |ci| 



Ica/M that 



depend on the value of /?. (iv) adjust the value of /3 us- 
ing a bisection method to find the zero of |cm-i-iP and 
repeat steps (ii)-(iv) until a desired tolerance is reached. 
In Fig. [1] we plot the numerically determined solution 
to the recurrence relation Eq. (|2.19p for AI — 29. The 
figure clearly illustrates the predicted symmetry about 
n = iM + l)/2. 

In Fig. [2] we plot the particle entanglement, 
Ep{\C^^^^^)), of the combined system plus ancilla 
against the total number of particles, M, in the ancilla. 



FIG. 2: The amount of particle entanglement Ep versus the 
number of particles in the reference ancilla, M, for a system 
consisting of a single shared particle. The ancilla is in the 
optimal reference state discussed in the text. 



The fact that Ep{\C^'^^+'^'))) > whereas Ep{\^^^^)) = 
i?p(|$(^))) = shows that the ancilla partially shields 
the entanglement of the system from the local U-SSR. 



D. An ansatz for the optimal reference 

The numerical solution depicted in Fig.[l]suggests that 
the analytical solution may be of a trigonometric form. 
To check, we now consider the ansatz for the solution to 
the recurrence relation of the form 



where 
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Here A, B, ^ and 6 are constant for a given value of M, 
i.e. a given number of particles in the ancilla. We note 
that the numerical solution in Fig. ([1]) comprises almost 
one full cycle of the cosine function. This implies ^ can 
be expressed conveniently in the form 



e = 



27r 



M + 2e 



(2.24) 



where e is an adjustable parameter. Moreover the sym- 
metry condition Eq. (jASP implies a condition on the value 
of 9 as follows. Replacing the left and right sides of 
Eq. (|2.2ip using the ansatz Eq. (|2.22p and the expres- 
sion for f in Eq. ()2.24p shows 
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which is satisfied by 



M + 2e 
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Hence the final form of our ansatz, given by Eqs. (|2.22p - 
(P:^ is 
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(2.27) 



and depends on just tliree parameters A, B and e. 

We now check the vahdity of the ansatz by substituting 
it into the recurrence relation Eq. (|2.19p . To simplify 
matters we use the form given by Eqs. (|2.22p and (|2.23p . 
The result of the substitution is 

(/3-l)|A-cos[z(n)]| 

= JA - cos[z(n + 1)] \{a~ cos[z{n - 1)] | 

+ |a - cos[2(7i)]| 

X |2A - cos[z(n + 1)] - cos[z{n - 1)] | , 

(2.28) 

which is independent of the normalization constant B. 
Using trigonometric identities and grouping terms it is 
easy to show that the right hand side of this equation 
can be expressed as 

RHS = iA^ + cos2(^) - 1 + [1 + 2 cos(0] 

x{ -2Acos[z(?i)]+cos^[z(n)]} . (2.29) 

Similarly the left hand side of Eq. ([2?22| is 

LHS = (/3-1) 1^2-271 cos[z(n)]+cos2[z(n)] I . (2.30) 

Comparing the coefhcients of cos[z(n)] and cos^[z(n)] in 
Eqs. (1^:^ and 1^:^ shows that the value of /3 that 
satisfies the recurrence relation is given by 



/3-l = 2cos(^) + l . 



Comparing the constant terms in these equations simi- 
larly shows that (/? - l)A^ = iA^ ~ iAv?{C)- These last 
two results imply that 
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1 + cos(C) 
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Eliminating cos(^) and sin(^) in Eq. (|2.32[) using 
Eq. (P3T|) then yields 



fi^AA^ 



(2.33) 



These results show that the ansatz does indeed satisfy the 
recurrence relation for particular values of the parameters 
A, B and e. 

To determine the actual values of the parameters we 
need to apply additional conditions. We first consider the 
boundary conditions |c_ip = and |cm+iP = 0. Setting 
n ~ —1 in the detailed form of the ansatz Eq. (|2.27p 
yields 
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(2.34) 



which gives an expression for A in terms of e as 

"27r(e-l)" 



A = cos 



M + 2e 



(2.35) 



The same result is obtained for n = M -\- 1 due to the 
symmetry condition Eq. (|2.2ip . Finally we note that the 
analytic solution |cip = |cop(/3 — 1) given in Eq. (|2.22p 
yields another expression involving A and e as follows. 
Using the ansatz Eq. (I2.27P to replace the coefficients 
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and |cip and Eq. ()2.33|) to replace (3 gives 
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A — cos(- 
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M + 2e' 

(2.36) 
This analysis shows that our ansatz, which depends on 
the values of just three parameters A, B and e, satisfies 
all the conditions for the optimal reference state provided 
Eqs. (|2.35p and (|2.36p are satisfied. Indeed the values 
of the parameters e and A are determined by solving 
Eqs. (|2.35p and (|2.36p simultaneously, and the value of 
B determined from the normalization condition Eq. (|2.7p . 
Once these values are determined, our ansatz Eq. (|2.27p 
provides an exact analytical expression for the optimal 
reference state for any value of M . 

We can also find approximate values for A, B and e in 
the large M regime. Substituting for A using Eq. (|2.35p 
and simplifying yields 



27re 



M- 



sm 



2e, 

27r(e 



27r 



M 



2e 



M + 2e 



4cos^ 



M + 2e 
27r(e 



1) 



M + 2e 



(2.31) which, in the large M regime, gives 



6 - - + o{Ar 



(2.37) 

Similarly, we find from Eq. (|2.35p in the same regime that 

A=l + OiM-^). (2.38) 

The normalization constant B is found by evaluating 



E. 



1 using the ansatz Eq. (|2.27p . i.e. 



M 



^=j: ^ 



n=0 



COS 



27r(n + e) 
M + 2e 



(2.39) 



For the case where M is odd we find, using the symmetry 
condition Eq. ^^, that 



B = 2J2 fin) 



where 



f{x) = < A — cos 



27r(x + e) 



M + 2e 



(2.40) 



(2.41) 



is a monotonically increasing function over the range 
from X = to X = ^'^t^ ■ We note that 



2 / dxf{x) >B>2 dxj{x - 1) , (2.42) 

Jo Jo 

that is, 

A{M + 1) + 1 + e > B > A{M + 1) - 3 + e , (2.43) 



and so 



B- 



1 



A{M + 1) 



+ 0{M-^) . 



(2.44) 



The same result is also found for the case where M is 
even. Hence in the large M regime the analytical form 
of the optimal reference state is 



1 



M + 1 

2 
M + 1 



1 



27r(n + 3/2) 



Af + 3 

7r(n + 3/2) 
Af + 3 



+ 0(Af-2) 
0{M-^) . (2.45) 



For comparison, in Fig. [T] we have also plotted as the 
solid curve the values of the probabilities |c„p given by 
the ansatz in Eq. (|2.27|) for M = 29 for the values e = | 
and A = 1, which are correct to 0{M^^). The figure 
shows that the ansatz with these approximate values of 
e and A is in relatively close agreement with the exact 
numerical solution (crosses) at this value of M. 



E. General nature of the optimal reference states 

To see why the optimal reference state has the gen- 
eral form shown in Fig. [1] and Eq. (|2.45p consider two 
reference states in opposing limiting cases. First, in the 
limit of a narrow particle number distribution only one 
coefficient is non-zero; it is straightforward to show using 
Eq. ([2?T6l) that the particle entanglement Ep{\C'^^'^+^'>)) 
is zero for this reference state. In the other extreme, 
the particle number distribution is uniform with |c„p = 
1/{M + 1) for n = 0, 1, . . . , A/; the right hand side of 
Eq. d^mi) is 1 - M/{M + 1) for this reference state. 
Clearly a reference system with broad particle num- 
ber distribution has an advantage for making parti- 
cle entanglement accessible. However, the reason for 
£'p(|C'^^^+^^)) being less than 1 bit in the broad particle 
number distribution case can be traced to the bound- 
ary conditions Eq. (|2.12p . This suggests that a distribu- 
tion which vanishes in near the boundaries n = and 
n = M also has an advantage. Evidently the particle 
number distribution of the optimal reference state bal- 
ances these two opposing effects, breadth of the distribu- 
tion with vanishing values near the boundaries, to maxi- 
mize ^p(|C(^^+i))). 

Our derivation of the optimal state fixes only the mod- 
ulus of the coefhcients c„ in Eq. (2.6). The arguments of 



the complex numbers, c„, are completely arbitrary. This 
is somewhat surprising given the role of the ancilla is pre- 
sumed to act as a phase reference. We now explore how 
the ability to act as a reference for the U-SSR differs from 
the ability to act as a phase reference. First note that 
the asymmetry of a state depenjis solely on the modulus 
of number state coefficients c„ 5]. For example, consider 
the state of the reference ancilla representing a uniform 
sharing of particles with arbitrary complex arguments 
given by 



|$(A'^)) 



AB 



^/M 



M 

f + 1,1^0 



'|n,Af -n) 



AB 



(2.46) 



where {6'„ : < ^n < 27r} are an unordered set of 
phase angles. Under the local group operation Tj^ik/S) = 
^-iNAkA ^ ^^ where A = 2tt/{N + 1) this state trans- 
forms to 



VM + 1 



E< 

n=0 



■'(<^'^-"M\n,M-n)AB 



(2.47) 



which is orthogonal to |$'-^^'')ab for integer k satisfying 
1 < k < N. This shows that |<i>^*'^'')AB is asymmetric 
with respect to U(l) and thus breaks the local U(l) sym- 
metry [39] . As this result is independent of the values of 
the phase angles {On} it confirms that only the m,oduli 
of the number state coefficients are important in terms of 
the U(1)-SSR. Moreover, Eq. ([^1^ shows that the par- 
ticle entanglement of a system consisting of one shared 
particle and an ancilla in reference state |$'-^^-')ab is in- 
dependent of the phases, {On}, of the reference state. 

Next consider the phase properties of the state 
|$^*^-')ab for various choices of the set of phase angles 
{On{- For this we use the Pegg-Barnett pha se formalism 
for physical states in the infinite-s limit [2J, [23, l2g . The 
joint phase probability density P$(0a,^b) for phase an- 
gles 6'a and 0^ which describe the phase operators (J)a 
and (/)B of the ancilla spatial modes at sites A and B, 
respectively, is given by 



PiOA^O^) 



{2^Y 



^e^("^*+'"«-)(n,TO|*(*^))AB 



(2.48) 
For the reference state |<i>'^*^))AB this can be rewritten as 



P(0A, ^b) = ^V^{Oa - Ob) (2.49) 



where 

Vii0A-OB) 



2-k{M + 1) 



M 

E' 

n=0 



,i[n(0A-eB)+6l„] 



(2.50) 



is the probability density for the phase difference (/)a — 4'b 
and the factor l/27r in Eq. (|2.48p is the uniform phase 




FIG. 3: The distribution of the relative phase difference for a 
random distribution of {9„}. 
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FIG. 4: The distribution of the relative phase difference where 
{9n} is linear in n. 



probability density for the phase operator of either spa- 
tial mode. The highly correlated nature of the phase 
difference stems from the sharing of a fixed number M 
of particles between the two modes. Figure [3] shows this 
phase distribution for the case where 9n = Trn whereas 
figure U] shows the phase distribution where 9n is random. 
Note that both phase distributions equally alleviate the 
U-SSR. Clearly there is a fundamental difference between 
the ability to act as a phase reference and the ability to 
act as a reference for the U-SSR, viz. the complex phase 
of the number state coefficients are crucial for the former 
but are unimportant for the latter. While this conclusion 
applies to the specific case of a single-shared particle sys- 
tem as considered in Eq. (|2.16p , we show later in Section 
IIVI that it also extends to arbitrary systems states. 



III. COMPARISON OF THE OPTIMAL 
REFERENCE FRAME TO OTHER STATES 

States with broad particle number distributions are 
known to have well-defined phase properties [2J, [2^, [2y, 
1271 . So it is interesting to compare our optimal reference 
states with classes of states that have been optimized 
for other phase dependent quantities. We define a figure 
of merit, D, based on the excess particle entanglement 
that is made accessible by the optimal reference state as 
follows: 



D = 



eMT'' 



-EA\C^ 



M+l)\ 



i?p[lC 



where 



|C(^+i)) = |^(i)),^|^(M)) 



(3.1) 



(3.2) 
(3.3) 



for single particle system state |\E''-^-'), optimal reference 
ancilla state |$opt ) ^^'^ arbitrary reference ancilla state 
\tp(^)'j. We now consider a number of states whose phase 
properties are important in some way. All the states are 
constructed on the Hilbert space Hr which is spanned by 
the number state basis {\n, M — n) : n = 0, 1, . . . , M}. 

The Berry- Wiseman phase optimized states of a two- 
mode optical field give the optimum phase shift estima- 
tion in a Mach-Zehnder apparatus for a fixed total num- 
ber of particles (photons) under ideal canonical phase 
measurements 281. These states have the form 



ex > sinM- — -!-\\n,M- 



n=0 



MH-2e 



(3.4) 



where M is the number of particles and e — \. The 
Summy-Pegg phase optimized states of a single mode op- 
tical field give the minimum phase variance for a fixed up- 
per bound in the particle (photon) number distribution 
[27| . Their two- mode version on the Hilbert space Hr is 
given by Eq. (|3.4p with a parameter value of e ~ 0.84 for 
M > 10. The variance of the optimized phase quantity 
for both classes of states scales as 1/M^. A class of states 
with less phase resolution is given by the single-mode co- 
herent states \a) oc ^„ a"|n)/vri! where \n) is the usual 
number state. These states have a phase variance which 
scales as l/|ap and are optimized in the sense that they 
app roximate number-phase minimum uncertainty states 
|25| for |a| 3> 1. To make a comparison with our opti- 
mal reference state we construct a two-mode version in 
the Hilbert space Hr with an analogous number state 
expansion to \a) as follows 




n, Af 



(3.5) 



The parameter that is an alogo us to the field amplitude 
a in \oi) has the value ^Mjl here. Thus in the large 
M regime the variance of the two-mode phase opera- 
tor scales as 2/M. We also consider reference states 
whose particle number probability distribution |c„p cor- 
responds to a binomial distribution, i.e. 



M 

E 



M 



p"(l_p)(M-«) 



1/2 



\n,M-n) (3.6) 



for which we set p = 0.5 to make the distribution sym- 
metric. These states are a two-mode generalization of 
Stoler ei aL's single mode binomial states [29]. Finally 
we include in our comparison the shared phase state, 



1 



y/M + 1 



M 

El 

n=0 



n, M 



(3.7) 



which is a two mode version of the sin gle- mode Pegg- 
Barnett phase state with zero phase [23, l2a] ■ This state 
belongs to the class of states on Hr with maximum asym- 
metry with respect to the particle number SSR Q. The 
change in the von Neumann entropy of these states due 
to the action of the SSR is the maximum value of log2 M 
and so these states play an important role in terms of 
breaking the symmetry represented by the SSR [5| . 

Fig. 3 shows that the figure of merit, D, of all states, 
at least after a certain value of M, are monotonically de- 
creasing with M. However, the decrease in D for increas- 
ing M for the Berry- Wiseman and the two-mode Summy- 
Pegg minimized phase variance states is far greater than 
the two-mode versions of the coherent and binomial 
states. The reason for this can be traced to our ear- 
lier observation in Sec. Ill El that broad particle-number 
distributions are an advantage in making the entangle- 
ment accessible. In particular the former two states have 
a particle-number standard deviation that scales as M 
whereas the two-mode coherent states have a particle- 
number standard deviation that scales as \/M. Com- 
pared to the maximum possible width M of the distri- 
bution, the former are relatively broad and the latter is 
relatively narrow. We note that the figure of merit D 
for the coherent and binomial states become closer as 
M increases as expected from the fact that the associ- 
ated particle-number distributions approach each other 
as M ^ oo. 



IV. OPTIMAL REFERENCE STATES FOR 
GENERAL PURE STATES 

We now generalize our results to the case where there 
are N particles in the system and, as before, M particles 
in the ancilla. Consider the tensor product of two arbi- 
trary states of the system and ancilla for this case given 

by 



|C( 



M+N}^ 



AB 



|*(^))ab ® |<i>(*'') 




FIG. 5: A comparison of the figure of merit D for the shared 
phase states (o), two- mode versions of the coherent states 
(*), binomial states (-I-) and Summy-Pegg phase optimized 
states (x), and Berry- Wiseman phase optimized states (D). 
D represents the relative effectiveness of a given reference 
state compared to the optimal reference state. The effective- 
ness of all states increases as the total number of particles M 
increases for large M. 



where the system and ancilla states are 

N 

|*('^))AB = ^d„|n,iV-n)AB , (4.1) 

n=0 



M 



|$^*'^)AB = ^C„|7Z,M-n)AB 



(4.2) 



n=0 



respectively. The coefficients c„ and d„ are subject to 
the normalization conditions 



N 



M 



^|d„p=^|c„p = l. 



(4.3) 



n=0 



m=0 



Neither the system nor the ancilla contain entanglement 
of particles i.e. £;p(|^(*='))ab) = £:p(|$(*'^))ab) = 0. To 
simplify the notation we shall omit the subscripts A and 
B from here onwards. The state of the combined system 
can be expressed as 

N M 

^(JM+N^^ = ^ ^ d„c„|n, N-n)® |m, M - m) 

n— ?72— 

N+M 

= E VSI^f ■'"^> (4.4) 



k=0 



where \C, 



(M+N)\_ 



is a state containing exactly k particles 



AB 



at site A and {N + M — k) at site B. Also p„ is a nor- 



10 



malization constant, i.e. 



I a 



{M+N)^ 



Pk 



N M 
V ^ n— m—0 

x\n,N - n) i^\m,M -m) (4.5) 

N M 
^y^ y^ \dnCra\'^6n+ra,k ■ (4.6) 

n— 771—0 



The calculation of Ep in Section [TTl for a system compris- 
ing one particle is easily extended to the case N particles 
as follows. We find that 



N+M 



£;p(|C(^-^+^))) = Y. PkEuiK 



(M+N) 



» 



N+M 



k=0 

N M 



fe=0 71=0 m=0 

N+M N M 






S,i+,7i,k log2 



\(ljiCn 



Pk Pk 

l^n^m I 



^ ^ ^ \dnCm\^Sn+m.k'^Og 



fc=0 n=0 771=0 



Pk 



(4.7) 



We again find that the particle entanglement is indepen- 
dent of the complex phases of the number state coef- 
ficients c„ of the reference state. This shows that our 
conclusion in Section III El for the special case of a single 
shared-particle system state, namely that a phase ref- 
erence and a reference for the U-SSR are fundamentally 
different, applies also to the general case of arbitrary sys- 
tem states. 



A. The infinite-particle ancilla state 

One may suspect that in the limit of M ^ oo, the an- 
cilla could resemble a classical reference and completely 
shield the system from the local particle number SSR. 
We now show this using the uniformly distributed an- 
cilla state with c„ = l/\/M + 1 for M ■> N . The state 
with exactly k particles at site A for A^ < fc < M is given 
by 



\C, 



and so 



(M+JV), 



1 ^ 

--=Ydn\n,N -n)® \{k-n),M - [k - n)) 



n^O 



i?p(|cf+"^))=i?M(|cf^+^^)), (4.8) 



and the probability "pk in Eq. (|4.6|) of finding this state is 
simply 



Pk 



N 

E 

71 = 



\dn 



1 



Af -f 1 M^\ 



(4.9) 



The contribution to £'p(|C^*^+^'')) from these states is 

M 



A^5:p.i^M(icr+"')) 



k=N 



M+l-N 
A/ + 1 



Em{\^^''^)) 



It is straightforward to show that the remaining 2N 
states |cf ^'^^^) for < fc < TV and M < fc < A^ -I- M 

have Euiici"'^^^)) < EmU^^''^)) and pk < 1/(M + 1) 
and so their contribution to Ep{\C^^~^^^)) is less than 



Y 



2N 
M + 1 



Euil^^''^)) 



These last two results give bounds on £'p(|C'^*^+^^)) as 

A<£;p(|C)) <X + Y. 
In the limit Af -> oo, A -> £;m(|*'^^)) and F ^ and 



so 



Ep{\C^"^+^^)) ^ Euil^^"^)) . 

The infinite-particle ancilla effectively shields the system 
from the particle number SSR and makes all of its en- 
tanglement accessible. We note that the class of perfect 
reference ancillae for the SSR include states that are poor 
phase references. 



B. Conditions for general optimal ancilla states 

We now use the Lagrange multiplier method to 
find conditions for the coefficients c„ of the opti- 
mal ancilla state for an arbitrary, but fixed, value of 
M. We want to maximize the particle entanglement 
i<;p(|C'(*^+^))) over |c„p subject to the normalization 



constraint ^ 



M 



1. The optimal ancilla state for 



any system state is given by the extremum of the auxil- 
iary function 



M 



F = i?p(|C(*^+^))) - ai{J2 IctiP - 1) (4.10) 



n=0 



where Ep{\C'-^''+'^^)) is given by Eq. ([iJ]) and ai is the 
Lagrange multiplier. Setting the derivative dF/d\c„ 



to 



zero we find with a little effort that 



o-i 



N 

E[-Ki 



\0g2{\dnCjn\^) + | ci„ |^ loga Pti-I-ti 



71=0 



(4.11) 
for ?n = 0, 1, . . . Af and where pk is given by Eq. (|4.6p . 

It is easy to show that the symmetry condition 
Eq. (|2.2ip does not hold in general. Consider, for ex- 
ample, the simplest case given by A = Af = 1, where 
both the system and ancilla contain a single shared par- 
ticle. Substituting these values into Eq. (14. lip yields for 
m — Q 



ai = |ciplog2 



|cidop-t- Icprfil 
Icodop 
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and for m 



ai = |coplog2 



|cido| 



IcodiP 



|codi|2 



where we have made use of the normaUzation Eq. (|4.3p to 
include an extra factor in the first line. Using this result 
we find the relationship between the coefficients |cmp for 
sequential values of the index m is 



The simultaneous solution of these equations is given by 

2 



1 



X 



/3|c„i+i 



|2(Af+l) _ 



p\c 



|2(M- 



1 - V™ 

-1) ^ ^j=0 



Em I 
7=0 |Ci 



where x = |ci/cop. This implies that the optimum refer- 
ence state for a non-symmetric system state, |co| 7^ |ci|, 
is also non-symmetric with |do| ^ \di\- 

The analytical treatment of system and ancilla states 
with arbitrary values of N and M, respectively, is beyond 
the scope the current paper. Instead we focus here on 
the situation where the system state is of a simple form, 
namely, a shared phase state with rf„ = l/y/N + 1, i.e. 
we set 



I3\ 



Cm — 1 



2(Af+l) _ 



= /3\cm\'^''+'^ 



i)l-S 



j=m I'-il 



Z^7=m I'-il 



(4.18) 



(4.19) 



We also find by setting specific values of m in Eq. (|4.17p 
that 



Ivi/W) 



^/N 



1 ^ 

^y\n,N-n) 



(4.12) 



We further restrain the problem to the case where N = 
M to make the problem tractable. 

A shared phase state has a symmetry between the sites 
in the sense that it is invariant under an interchange or 
site labels, i.e. AB h^ BA. The corresponding optimal 
reference state has the same symmetry with coefficients 
of the form 



M n 

M Af-1 M 



(4.20) 



(4.21) 



j=o 



n— j^n 



and so |co| 



CM 



Using this result together with 



\Cm I 



\CM-r 



(4.13) 



for TO = 0, 1, ... , M. The proof of this is as follows. We 
find from Eqs. ()4.1ip and (|4.6p the set of equations 



Eqs. (|4.18p and (|4.19p with the specific values of m = 
and m — M, respectively, shows that |cip — |c7\/_ip. 
Repeating this analysis for sequential values of m com- 
pletes the proof of the symmetry in Eq. (|4.13p . 

Another property of the optimal state is that the dis- 
tribution |cmP is unimodal which can be seen as follows. 
First we define the difference Am+i,m as 



^771-1-1,771 l^in+l 



|2(M+1) _ I |2(M+1) 



(4.22) 



ai 



log2 



for TO = 0, 1, . . . , M where 

1 

Pk 



N- 



N 
ri=0 



log^ Pn+m (4.14) which, from Eq. (|4.18p, is given by 



A^+i,„^|c^p(^'^+^)( ' '^T°L )■ (4-23) 



j:T=o\cA' 



N- 



I E 1^.1^ 



(4.15) 



j = k-N 



Here, for compactness, we have extended the set of coef- 
ficients c„ with the values c„ = for n < and n > M. 
Taking the exponential of both sides of Eq. (I4.14p yields 



From the symmetry property Eq. (|4.13p of the optimal 
state and the normalization condition Eq. (j4.3p we know 
that 



[M/2] 

E 1^" 

m=0 



M 

E 

= [Af/2]H 



(4.24) 



A'' n-\-m 



/3iciiir+^) = n E 



n— j—n+m—N 



(4.16) 



Omitting the zero terms |cjp for j < and j > M and 
setting N = M we find 



where [M/2] is the largest integer satisfying [Af/2] < 
M/2. Next we use this result to determine the posi- 
tivity of the differences A,„+i .,„. Consider if the ex- 
pression 1 — 2^"1q |cjP in Eq. (|4.23p were zero for a 
value of 771 less than [A// 2]. This would imply that 



M — 771 71 + 771 AI AI 

; i^.f) n ( E 1^.'^ 

n—O j — n—AI — m j—n+m — AI 

M n m M 



"m+l| , 



|Cm+2h 



P\cm\'^^-''^=U (Ei^.i^) n ( E 

n=0 j=0 
M n 

n (E 1^.1^) n(E 1^.1% (4-17) 



-[Af/2] I 



and by Eq. KTS^ that 



^m+i.m ~ 0, which would mean that \cm\ = |cm+i| = 
and so 1 — 2^™T) \cj\^ = 0. It follows, by in- 



duction, that the coefficients 



would be zero for 



j = 0, 1, . . . , [M/2] which contradicts Eq. ((424)) . Hence 



71—771 j—0 71—0 j—n 



l-2E"o 



only for m — [M/2] and so from 
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FIG. 6: Comparison of the particle number distributions for 
the optimal reference state |"1>^ ^ ) (x) and the trial state 

I'^lriS) (— ) for Af = M = 29. The comparison shows that 
the optimal reference state is very well approximated by the 
trial state. Moreover, a comparison with the optimal state 
for A'^ = 1,M = 29(— — ) shows that the optimal particle 
number distribution broadens out as the number of particles 
in the system, A'^, increases. Finally we note that distribution 
for the optimal reference state ( x ) satisfies the symmetry and 
unimodal properties derived analytically. 



Eq. (li:^ Aj+i^j > for j = 0, 1, . . . , {[M/2] - 1). Tak- 
ing account of the symmetry Eq. (|4.13|) then shows that 
the distribution |c„p is monotonically increasing over the 
index m for m =^ 0, 1, ... , [A//2] and monotonicaUy de- 
creasing for TO = [M/2], {[M/2] + l),...,M. This com- 
pletes the proof of the unimodal property. 

In summary, the optimal state for a system in the 
shared phase state Eq. (|4.12p and a reference ancilla with 
an equal maximum number of particles, i.e. with M — N, 
has coefficients |c,„| that are symmetric about a modal 
point and are monotonically increasing before this point 
and are monotonically decreasing after it. These proper- 
ties are reminiscent of the ansatz used for a single shared 
particle case which suggests that the optimal state would 
be well approximated by a trial state |(/'J:,.iai) which is of 
the form Eq. (|2.27[) for a suitable choice of the parame- 
ters A and e, and a range of M values. Exact analytical 
solutions or better approximations are beyond the scope 
of this paper. However, one can solve the N — M case 
numerically using either Eq. (|4.18p or Eq. (|4.19p . Fig. [g] 
shows the particle number distribution for the exact op- 
timal reference state with N — M — 29 compared with 
that of the trial state for the parameter values A — 1.9 
and e = 8.9. The figure confirms our intuition that the 
exact solution in the N = M case can be approximated 
by a state which is of the form Eq. (|2.27p . The inner 

product ('&opt I'/'triai) "^^ *^^ exact solution with the trial 



state differs from one by 6 x 10 



V. DISCUSSION 

We have focused on the situation where the ancilla 
state |$(*^)) has a particularly simple form, as given by 
Eq. (|2.6p and Eq. (|4.2p . corresponding to a single spa- 
tial mode at each site. Such states automatically satisfy 
the requirement EpQ^'^^^)) = 0, i.e. that the ancilla 
does not contain particle entanglement. Nevertheless our 
analysis generalizes quite easily to situations with multi- 
ple spatial modes and multiple ancillae. The most gen- 
eral pure ancilla state of M particles which contains no 
particle entanglement has the form 



1$ 



(A/)\ 



M 

E 



c™|(/)(to))a ® [(p{M - to))i 



(5.1) 



where |x("-))a" ^^r x = or 1^9, is an arbitrary normal- 
ized state of exactly n particles at site /^ in the second 
quantization formalism. For example, if the part of the 
ancilla at site /^ involves the 3 mutually orthogonal spa- 
tial modes |-)^ ® [)fi® [)in then the general expression 
for Ix("-))m is 

IxH)m = E fi-]'k[i)A ® |J)a ® [k)A5i+j+k,n (5.2) 

where fi,j,k are normalized coefficients. Notice that 
the set of states at each site are orthonormal, viz. 
m(x("-)Ix("^))a' = ^n,m- To incorporate these more gen- 
eral ancilla states in our previous analysis we need only 
replace the ancilla states [n)^ in the above with |(/)^(n))^. 
However, as our analysis used only the orthogonality 
property of the set of |n)^ and not its detailed struc- 
ture, the coefhcients |cmP for the optimal ancilla state 
would be unchanged. Hence the solutions derived above 
represent classes of optimal ancilla states for arbitrary 
choices of coefficients of the kind fi.j^k m Eq. (|5.2p . 

Moreover consider an ancilla which comprises compo- 
nent subsystems of the form 



where M = Mi + M2. The 
£;p(|$'(*'^))) = implies that 1$'^*^^; 
in the form 



(5.3) 

requirement that 
can be expressed 



w 



(A/), 



M 
m=0 



\cj){m))A®[v{M~m))i 



(5.4) 



where the tensor product form of Eq. (|5.3p restricts the 
variability of the coefficients cj„. This is a further con- 
straint on the coefficients, in addition to the normaliza- 
tion condition, and so the entanglement made accessible 
by a composite ancilla cannot exceed that given by a 
single component ancilla composed of a maximum of M 
particles. Hence the use of multiple ancillae in a tensor 
product state does not give any advantage in optimizing 
the accessible entanglement. 
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VI. CONCLUSION 

Quantum reference frames are of particular theoretical 
interest and the phenomena of unlocking entanglement 
has proven to be a very interesting application. Our ana- 
lytical techniques have proven to be fruitful in producing 
exact solutions in the case of a single shared particle and 
guidelines for sensible approximations for system states 
of large N. Much work is needed to produce exact ana- 
lytical solutions for more general system states. However, 
the results extracted seem to establish the intuition that 
pure states with a broad particle-number distribution act 
as good reference frames for the particle number SSR. In- 
deed, the optimal reference state in the infinite-particle 
limit M — > oo makes all the entanglement in the single 
shared particle accessible. 
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APPENDIX A: SOLUTION AS A SERIES OF 
POLYNOMIALS IN /? 

Here we explore the nature of the iterative solution 
of the recurrence relation Eq. (|2.20p for the case where 
the system consists of a single particle which is equally 
shared between the sites. Setting n = and using the 
boundary condition |c_ip = gives an expression for 
|cip in terms of |cop. Next, by setting n = 1 we find 
an expression for |cip and so on. Continuing in this way 
yields a solution to Eq. (|2.20p that is expressed in terms 
of just f3 and |coP as 



|c„P=P„(/3)|co|^ 



(Al) 



where P„ (/?) are polynomials in /3 of order n given by 



Po{0] 


= 1 


pm 


= /3-l 


P2{P) 


= (/5 - 1) (/5 - 2) 


pm 


= (/32-3/3-fl)(/3-2) 


pm 


= (/3-3)(/3-l)(/32-3/3+l) 


pm) 


= (/3^-5/32-h6/3-l)(/3-3)(/3-l) 


pm) 


= (/33 - 6/3^ + 10^ - 4) (/33 - 5/3^ + 6/3- 1) 




(A2) 



etc. Due to the normalization condition Eq. (|2.7p . we 
have the property that 



M 



E^»(/3) 



n=0 



|co| 



(A3) 



Alternatively, rearranging the recurrence relation 
Eq. (jTTg)) as 



|c„-i| 



|c„P-H|c„+iP 



(A4) 



leads to a different form of the solution. For this we 
set n — M and use the boundary condition c^vz+i = 
with Eq. (|A4[) to get an expression for |ca/_iP in terms 
of |cA/p. Continuing in this way with reducing values 
of n we find that any coefficient |cM-nP can be written 
in terms of polynomials in P and | cm \ ^ • The symmetry 
of Eq. (I2.19P with respect to interchanging n with m = 
M — n ensures that the polynomials are the same as those 
appearing in Eq. (jAip . i.e. 

\cM-nf = PnmCAlf ■ (A5) 

Substituting n = M into Eq. ([A5|l and (|AT1) yields 

|coP = Pm(/3)|cmP , |CA/P = Pa/(/3)|coP (A6) 

respectively. The simultaneous solution to Eqs. (jA6p is 
Pm(/3) = 1, |coP = |ca/P (A7) 

and so from Eqs. (|Aip and (lASP we find that 

|c„P=|cA/-nP (A8) 

for n — 0,1,..., M. Thus the solution which represents 
the optimal reference state is symmetric in this sense. 



APPENDIX B: PHASE PROPERTIES OF THE 
STATE |$i*^')ab 

The unitary operator Ua = e'^*^^*"^''' in the defini- 
tion of the state \^[ ' ')ab in Eq. (??) represents a Kerr- 
like nonlinear interaction. The action of this operator on 
coherent states has been discussed extensively in the lit- 
erature (see e.g. Refs. [30,[311[32,[3SlM])- However, in 
our case the operator {/a acts on the shared phase state 
1^0 1)ab in Eq. (??) and so requires a minor extension 
of previous work. For this we use a previous result to 
derive a form of the operator JJ^ that allows us to infer 
the pertinent phase properties of \^\ )ab- 

One of us has previously shown that 34 1 



A'-l 



> Cfeexp(m0fc) — exp[— M9n(n — 1)] (Bl) 



fe=0 



1^' 



— y ^ exp[— z'i?n(n — 1)] exp(m^j) (B2) 



K ^^ 

n=0 



where 



K 



'>k =9Vi 



27rfc 



mod2} 



K 



(B3) 

(B4) 
(B5) 
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and J and K are positive integers sharing no common 
factors. Expanding [/ in terms of the Fock basis then 
yields 

oo 

f7A = ^exp[-iM"-l)]|n)AA("| (B6) 

oo K-1 

= ^ ^ Cfcexp(m(/)fc)|r^)AA(?^| (B7) 

K-1 

= ^ Cfc exp(iiVA</)fc) (B8) 

fc=0 

which is in the form of a series of phase shifting operators 



exp{iN(j)k)- The case discussed in the text is for the 
choice J = 1 and X = 2 for which -d = 7r/2, 0o = 7r/2, 
4>i — 3tt/2, Co = (1 + J)/2, and ci = (1 — i)/2, and so 



M 



l + i 



;^' Ve— /^|n,M-n)AB 

^ n— 

+ ^717^E/"^"^^'I"'^-^)ab .(B9) 



n=0 



Using Eq. ()2.50|) then shows that Eq. (??) is the phase 
difference probabihty density corresponding to this state. 
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